A Non-Oracle Quantum Search Algorithm and Its Experimental Implementation 
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Grover's algorithm has achieved great success. But quantum search algorithms still are not 
complete algorithms because of Grover's Oracle. We concerned on this problem and present a new 
quantum search algorithm in adiabatic model without Oracle. We analyze the general difBculties 
in quantum search algorithms and show how to solve them in the present algorithm. As well this 
algorithm could deal with both single-solution and multi-solution searches without modification. 
We also implement this algorithm on NMR quantum computer. It is the first experiment which 
perform a real quantum database search rather than a marked-state search. 



PACS numbers: 03.67.Lx, 89.70.-a, 03.65.-w 

Quantum computation is a promising way to solve clas- 
sical hard problems. Although large-scale quantum hard- 
ware has yet been built, quantum computation model in 
analog to classical circuit is well developed during the 
last few years. Based on this model, several quantum 
algorithms have been designed to perform classical algo- 
rithms with remarkable specdups. The most splendid one 
among these is Shor's Algorithm[l], which can factorizc 
a big number using a running time only polynomial in 
the size of the number, while all known classical algo- 
rithms need a exponential time[2]. Another important 
algorithm[3], named after its inventor Grover, concerns 
the problem of searching for a required item in a un- 
sorted database. One common example for this unsorted 
database search is to find a person's name in a phone 
book (the items are sorted by names) with only knowing 
his phone number. Classically, the only way to achieve 
this is brute-force search[4] which for N entries in the 
phone book requires an average of ^ quires. However, if 
the information is stored in a quantum database, to find 
the right name with Grover's algorithm costs only a time 
of order -y/iV, providing a quadratic speedup. 

While quantum algorithms arc presented in standard 
circuit model(i.e., using a sequence of discrete quantum 
gates), a new model of quantum computation show up 
where the states of quantum computer evolves continu- 
ously and adiabatically under a certain time-dependent 
Hamiltonian. This new idea was firstly brought out 
by Far hi and co-workers [5]. In this new computation 
model, a problem Hamiltonian is well designed whose 
ground state encodes the unknown solution to the prob- 
lem. Then this adiabatic evolution can be used to switch 
gradually from an initial Hamiltonian whose ground state 
is known, to the problem Hamiltonian. If this evolution 
evolves slowly enough, the system will stay near its in- 
stantaneous ground state[G]. So in the end of evolution, 
the system will on the solution state of the problem. This 
method has been applied to the database search prob- 
lem. However, this adiabatic search algorithm results in 



a complexity of order N, which is the same order with 
classical algorithms. More recently, Roland and Cerf[7] 
improved the performance of adiabatic search to order 
a/TV, the same with Grover's algorithm, by applying adi- 
abatic evolution locally. 

Although these quantum search algorithms seems bril- 
liant as they have already done, they are still incom- 
plete algorithms. Grover's algorithm utilized a Oracle 
(i.e., a blackbox) , which gets an input state |i), checks 
the quantum database, and changes the state to — |i) if 
the i-th value in the database satisfies the search condi- 
tion and does nothing otherwise. It is easy to implement 
such operations in classical database cases, but up to now 
there's no efficient universal method to design this Ora- 
cle in quantum circuit. And in adiabatic algorithms, the 
solution of the problem is encoded to the problem Hamil- 
tonian. Since the mechanics of the Oracle remains un- 
known, the encoding process of the Hamiltonian in the 
adiabatic algorithm is unclear. Instead, just like what 
previous experiments [8, 9, 10, 11] of Grover's algorithm 
did, the adiabatic search algorithm forms the Hamilto- 
nian directly from the solution state, which means we 
have to know the state in prior and then perform a al- 
gorithm to show it. Obviously this marked-state search 
algorithm is not a real database search. Thus the main 
problem with current quantum search algorithms is the 
existence of Grover's Oracle. 

In this article, we present a new adiabatic algorithm for 
quantum search. By encoding the database to quantum 
format and forming the problem Hamiltonian from tar- 
get value, this adiabatic search algorithm solves Grover's 
problem without Oracles. Furthermore, we experimen- 
tally implement this non-Oracle quantum search algo- 
rithm in NMR quantum computer. Because of the rea- 
sons mentioned before, this is the first time implementing 
a real quantum unsorted database search in experiment. 
Wc also analyze the general difficulties in quantum search 
and show how to solve them in our algorithm. 

As a new quantum computation model, adiabatic algo- 
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rithm was brought out by Fahi et al. [5] and soon became 
a rapidly expanding field. This new computing model 
relies on the Adiabatic Theorem which states as follows: 

Adiabatic Theorem. A physical system remains in its 
instantaneous eigenstate if a given perturbation is acting 
on it slowly enough and if there is a gap between the 
eigenvalue and the rest of the Hamiltonian's spectrum. 

The idea of adiabatic quantum computation is 
straightforward. First find a complex Hamiltonian whose 
ground state describes the solution to the problem of in- 
terest, Next, prepare a system with a simple Hamilto- 
nian and initialized to the ground state. Finally, the 
simple Hamiltonian adiabatically switches to the com- 
plex Hamiltonian. According to adiabatic theorem, the 
system stays in the ground state, so in the end the state 
of the system describes the solution to the problem. The 
time dependent Hamiltonian is usually constructed as fol- 
lows, 

H{t) - [1 - s{t)]H, + s{t)Hp, (1) 

where Hi is the initial Hamiltonian whose ground state is 
easy to know and Hp is the complex problem Hamiltonian 
whose ground state describe the solution to our problem, 
and the monotonic function s{t) fulfills s(0) — and 
s{T) = 1. 

Several adiabatic algorithms have been designed for 
solving computational hard problems[12, 1.3]. And a sim- 
ple proof has been given to show that adiabatic model is 
equivalent to circuit mode in quantum computation[14]. 
Moreover, since adiabatic computation only involves the 
ground state, it keeps the system at a low temperature. 
Thus the system appears lower sensitive to some per- 
turbations and have a improved robustness against de- 
phasing, environmental noise and some unitary control 
errors[15, 16]. 

As mentioned before, the key part of an adiabatic al- 
gorithm is how to describe the solution to a specific prob- 
lem in the problem Hamiltonian. Here let's focus on the 
database search problem. To be simplified, we now con- 
sider a phone book which contains N (assume N ^ 2") 
entries with each entry a pair of telephone number and 
person's name. Usually, the entries are sorted by name. 
The database search problem here is to find a specific 
name in the book whose telephone number is given. To 
solve this problem in our model, the names are encoded 
to 71-qubit states and the phone numbers represented as 
integers (in fact, any real numbers are permit). An ex- 
ample for TV = 4 is shown in Table I. We encode the 
names and the phone numbers as in Table II. Thus the 
database could be stored by the state-integer pairs like 
{(|0),4),(|1),3),(|2),1),(|3),2)}. If we want to find the 
name which connect to the number 3601003 which is en- 
coded as 3, state |1) should be returned from our quan- 
tum search machine. 



TABLE I: phone book example 





Name 


Number 




Alex 


3601004 




Bob 


3601003 




Cherry 3601001 




David 


3601002 


TABLE II: encodor table 


Name 


State 


Number Integer 


Alex 


|0) 


3601001 1 


Bob 


|1) 


3601002 2 


Cherry 


|2) 


3601003 3 


David 


|3) 


3601004 4 



After encoding classical database to quantum 
database, next step is to design the problem Hamilto- 
nian Hp. For this problem. Hp is. 

N-l 

Hp ~ ( numberi\namei){namei\—target-I®"')'^ (2) 

i=Q 

where target is the code of the phone number which we 
want to search for. Obviously, the ground state of Hp 
is the state which connects to the target number. In 
general, for an encoded database where entries are pairs 
as ii^valuei) and sorted by i, we write Hp as 

JV-l 

Hp = valuei\i){i\- target ■ I'^'^f 

^ - 2 ■ target ■ V + target'^ , (3) 

where T) = ^f^Q^ valuei\i){i\ is the database operator 
and could be formulated separately. 

Next, we will choose an initial Hamiltonian Hi. 
Convcntionlly[12], Hi is chosen to be noncommutative 
with Hp to avoid crossing of energy levels. Thus we write 
Hi as: 

H, = ,9(^2 + + • + (4) 

which means the qubits coulpling with a magnetic field 
at the x-dircction and the coupling strength is g. the 
gound state of this Hamiltonian is simple and they are, 

|0(«-1)) - |l("-l)) |0("-2)\ _ |l(n-2)N 

10(0)) - |l(o)) 
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FIG. 1: Energy diagram for searching an item {i.e., 3601003) 
in the example database (see Table 1). During the adiabatic 
evolution, the state will remain on the lowest energy level. 



where b{j) is the Hamming distance between j and 0. 

In the adiabatic evohition, the system Hamiltonian in- 
terpolates from Hi to Hp {i.e., see Eq 1) and the state 
of the system evolves according to the Schrodinger equa- 
tion: 



H{t)\m) 



(6) 
(7) 



If this evolution acts slow enough {i.e., the total evolu- 
tion time T is large enough), the Adiabatic Theorem 
ensures the system will always stay on the ground state 
of H{t) and in the end the solution of our problem will 
show up. 

Again, we take the phone book in Table I as example. 
If we want to find number 3601002 in the database and 
using the encoder in Table II, we will get the Hamiltoni- 
ans as follows, 

H, = g{al+al) 



Hr, 



(4|0)(0|+3|l)(l| + l|2)(2| + 2|3)(3|-2/)2 



I' 



.(1) 



(8) 



And the eigenvalues of time dependent Hamiltonian H{t) 
(see EQ.(l) ) are plotted in FIG.l. By the adiabatic the- 
orem, the state will stay on the lowest energy level during 
the adiabatic evolution. And finally we will get the state 
on the basis |3). After measurement and decoding we will 
get the name connecting to the number 3601002 which is 
David. 

To demonstrate this Non-Oracle search algorithm, we 
selected ^'^C-labeled CHCI3 as a physical system for our 



experiments. The two qubits are represented by -'^^C 
and ^H. Its natural Hamiltonian in the multiply rotat- 
ing frame is 



(9) 



where lui and UJ2 are Larmour frequencies, J is the 
spin-spin coupling constant J = 2\A.bHz. Experi- 
ments were performed at room temperature using a stan- 
dard 400MHz NMR spectrometer (AV-400 Bruker instru- 
ment). 
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FIG. 2: The quantum network for adiabatic search using 
NMR interferometry. The input state is |00). 6s = 0.95 * 
(1 - s/10), Ts = 0.95 * (s/10)/7rJ. 

The experiments was devided into three parts, shown 
as Fig. 2: the first part consists of preparation of the state 
of the initial Hamiltonian. The second part is the adia- 
batic evolution, and the third part is the tomography of 
the resultant state. To prepare the initial ground state 
[Eq. 5], we first created a pseudopure state (PPS) [17, 18] 
Poo = i^I+e|00)(00|, where e « lO'^ describes the ther- 
mal polarization of the system and I is a unit matrix, 
using the method of spatial averaging. Then 7r/2 pulses 
along the — y axis was applied to prepare the ground 
state. 

Discretizing a continuous Hamiltonian is a straight- 
forward process and changes the run time T of the 
adiabatic algorithm only polynomially. Simply, let the 
discrete time Hamiltonian H{s) be a linear interpola- 
tion from some beginning Hamiltonian H{Q) = Hq to 
some final problem Hamiltonian H{S) ~ Hi such that 
H{s) = {s/S)Hi + {1- s/S)Ho. The unitary evolution 
of the discrete algorithm can be written as 



-iH{s)T 



(10) 



where t = T/{S + 1), T is the total duration of the 
adiabatic passage and 5 -I- 1 is the total number of dis- 
cretization steps. When both T, S ^ 00 and r ^ 0, the 
adiabatic limit is achieved. 

For our example shown in Eq.8, an optimized set of 
parameters was set as 5 = 1, T = 10.45 and S = 10, 
so r = T/{S + 1) = 0.95. This set of parameters yields 
an adiabatic evolution that finds the solution in a rela- 
tively efficient way. Using the Trotter formula, we can 
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approximate Us to second order 



g-»{l-f)i/oig- 



Us « U's 

(11) 

the fidelity of Us U^ is all above 0.996 and over- 
all fidelity is 0.991. For the implementation of U'^, 
e~'(i~f)^oi- can be simply realized using a 6s pulses 
around x axis for both H and C nuclei, 9s = r* (1 — s/ S), 
shown in Fig. 2. And the evolution under Hi can be sim- 
ulated by a free evolution under the Hamiltonian Tisys , 
the identity term of Hi docs not cause any evolution of 
the state and so it can be omitted, = r * {s/S)/ttJ. 
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FIG. 3: The tomography of theoretically expected and exper- 
imentally obtained density matrices for the search states in 
adiabatic search algorithm. The density matrices consist of 
just a real term on the diagonal corresponding to the popu- 
lation of the state that has been searched. 

The third stage of the experiment is the tomography 
of the final density matrix after the adiabatic evolution. 
The result was shown as Fig. 3. Theoretically, the four 
state |00), |01), |10) and |11) should be find at the proba- 
bility of 0, 0.014, 0.014 and 0.972. Our experiments show 
that the probability is 0.037, 0.032, 0.006 and 0.925. The 
fidelity of the experiment is 0.985. 

The errors in adiabatic algorithms may be caused by 
three parts. Firstly, the total time of evolution in adia- 
batic algorithms should be infinite. Actually the evolu- 
tion is terminated when the state is supposed to reach 
our expected high probability. Secondly, the error is due 
to neglect of C(t^) terms in the Trotter Formula (Eq. 
11). The third part of error is due to decoherence effects 
of the NMR system and imperfect pulses. 

Unlike previous standard quantum algorithms only us- 
ing qubits as registers to store information, our algorithm 
represents the value field by the strength of interactions 
in the operator and the index field by qubits. This is 
because that if both the fields are represented in qubits, 
2n qubits are needed for a database with N items, which 
result in the failure that the optimal running time would 
be scaled from order y/N to N[7], the same performance 
as classical algorithms. Since the construction function 
is simply quadratic, the interaction strength in problem 
Hamiltonian grows with the database's size. For further 
consideration, the algorithm may be improve to suppress 
the interaction strength in the problem Hamiltonian by 



choosing a better construction function in Eq. 3. 

All practical quantum search algorithms must face the 
problem that the database is unsorted, thus quantum 
operators would traverse all the items in the database to 
learn the complete information and after measurement 
all information in the states are destroyed. It is a hard 
problem to efficiently implement the quantum operators 
concerning the database. The first effort was reported 
by Ju and coworkers [4] when they tried to implement 
Grover's Oracle in quantum circuit. They have to spend 

steps to build the relation in the database to the states 
for each query, such that the circuit design is not efficient. 
In this algorithm, we describe the database information 
in a single operator(I? in Eq.3), thus this operator may 
be analyzed and formulate separately for each database. 
On the other side, if it could not be formulated efficiently 
in some cases, approximate implementation is another 
possible solution according to recent works on geometric 
quantum computation [19], of which detailed considera- 
tion is beyond the scope of this article. 

As an end of this section, we will give a simple anal- 
ysis of the multi-solution search case in our algorithm. 
If there're m > 1 entries in the database satisfying 
the search condition, the problem Hamiltonian will have 
ground states with m times degenerated. And because of 
symmetry, the state would finally evolute to an average 
superposition of all the ground states. Thus without any 
modification, our algorithm could also deal with multi- 
solution search. 

To be concluded, we present a new kind of adiabatic 
search algorithm to solve Grover's problem without Ora- 
cles and give a demonstrative experiment on NMR quan- 
tum computer. The result of experiment agrees well with 
theoretical expectation. This is a new style of quantum 
search algorithm which utilize both quantum registers 
and interaction strength to store information. This al- 
gorithm aims at general difficulties of quantum search 
algorithms and give a promising way to solve them uti- 
mately. 
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